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Direct numerical simulations of three-dimensional (3D) homogeneous turbulence under 
rapid rigid rotation are conducted to examine the predictions of resonant wave theory for 
both small Rossby number and large Reynolds number. The theory predicts that "slow 
modes" of the velocity, with zero wavenumber parallel to the rotation axis (k z = 0), 
will decouple from the remaining "fast modes" and solve an autonomous system of two- 
dimensional (2D) Navier-Stokes equations for the horizontal velocity components, normal 
to the rotation axis, and a 2D passive scalar equation for the vertical velocity component, 
parallel to the rotation axis. The Navier-Stokes equation for three-dimensional rotating 
turbulence is solved in a 128 3 mesh after being diagonalized via "helical decomposi- 
tion" into normal modes of the Coriolis term. A force supplies constant energy input 
at intermediate scales. To verify the theory, we set up parallel simulations for the 2D 
Navier-Stokes equation and 2D passive scalar equation to compare them with the slow- 
mode dynamics of the 3D rotating turbulence. The simulation results reveal that there 
is a clear inverse energy cascade to the large scales, as predicted by 2D Navier-Stokes 
equations for resonant interactions of slow modes. As the rotation rate increases, the 
vertically-averaged horizontal velocity field from 3D Navier-Stokes converges to the ve- 
locity field from 2D Navier-Stokes, as measured by the energy in their difference field. 
Likewise, the vertically-averaged vertical velocity from 3D Navier-Stokes converges to a 
solution of the 2D passive scalar equation. The slow-mode energy spectrum approaches 
k h 5 ^ 3 , where k% is the horizontal wavenumber, and energy flux becomes closer to con- 
stant, as in 2D, the greater the rotation rate. Furthermore, the energy flux directly into 
small wave numbers in the k z — plane from non-resonant interactions decreases, while 
fast-mode energy concentrates closer to that plane. The simulations are consistent with 
an increasingly dominant role of resonant triads for more rapid rotation. 



1. Introduction 

Large-scale flows in oceans and the atmosphere are greatly affected by the earth's 
rotation and are known to be quasi- two-dimensional. Rotation also plays an important 
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role in many engineering flows, e.g. high-Reynolds turbulent flows in turbo-machinery 
and rotating channels. When a fluid is under rotation, the Coriolis force is introduced into 
the momentum equation and competes with the nonlinear force. A dimensionless number, 
the Rossby number Ro, may be defined as the ratio of the magnitude of the incrtial term 
to the Coriolis force. In a rapidly rotating fluid, a mathematical limit Ro — > is taken 
of the three-dimensional Navier-Stokes equations. There is a general expectation that 
the fluid should become two-dimensional in this limit. For example, the classic Taylor- 
Proudman theorem implies that the Coriolis force will align vortex tubes parallel to the 
rotation axis in steady-state, rapidly-rotating flows (Greenspan (1968)). 

To take into account the effect of rapid rotation on dynamical evolution a resonant 
wave theory has been applied (Greenspan (1968), Greenspan (1969), Waleffe (1993)). 
Similar theories were first developed for gravity waves in geophysical fluid flows (Phillips 
(I960)) and have since been widely invoked elsewhere. For excellent reviews, see Phillips 
(1981) and Craik (1985). According to this approach, the fluid velocity may be regarded 
in the limit of small Rossby number as a superposition of inertial waves with a large 
characteristic frequency which are modulated on a longer, slow time-scale. An averaged 
equation is derived from weakly nonlinear theory for the slow-time motion. This equation 
explains enhanced energy transfer from small scales to large scales by the resonant tri- 
adic interaction of inertial waves (Greenspan (1968), Greenspan (1969), Waleffe (1993)). 
Using an "instability hypothesis", Waleffe (1993) argued further that resonant triadic 
interactions should drive the flow to become quasi-two-dimensional. More recently, this 
type of resonant wave theory has been put on a sounder mathematical footing. Embid 
& Majda (1996), Embid & Majda (1998) and Majda & Embid (1998) have derived sim- 
ilar "averaged equations" in a rigorous asymptotic limit over a fixed time interval for a 
general class of geophysical fluids problems with fast wave dynamics. For the particular 
case of rotating incompressible fluids, it has been shown that the averaged equations 
contain as a subset the 2D Navier-Stokes equations for the vertically-averaged velocity 
fields (Mahalov & Zhou (1996), Babin, Mahalov & Nicolacnko (1996)). 

However, in turbulent flow under rapid rotation not only is the Rossby number Ro 
small, but also the Reynolds number Re is large. As a consequence, eddy motions are 
excited on a wide range of length scales. At any given wavenumber k there are thus 
at least two distinct time scales, the rotation time scale rn ~ ^ and the nonlinear 
time scale T non (k) ~ (k 3 E(k))~ 1 ^ 2 (where E(k) is the energy spectrum). The validity of 
the wave resonance theories depend upon tq being shorter than all other time scales in 
the problem. However, for very large Reynolds numbers, there will be a range of high 
wavenumbers k where instead T non (k) <C tq- Thus, the resonant wave theory is likely 
to be valid for turbulent flows only very non-uniformly in wavenumber, if at all. The 
existing mathematical derivations of the theory (Embid & Majda (1996), Embid & Majda 
(1998), Majda & Embid (1998), Mahalov & Zhou (1996), Babin, Mahalov & Nicolaenko 
(1996)) are only usefully valid for low Reynolds numbers. Indeed, error fields for the 
resonant wave approximation are estimated in these proofs by "Sobolev norms" that 
get most of their contribution from high-wavenumbers. The errors are therefore shown, 
by present theorems, only to have upper bounds that grow rapidly with the Reynolds 
number Re. While the error bounds also decrease with Rossby number Ro, arguments 
based on comparing time-scales, like those above, suggest that resonant triads will be 
selected latest at the highest wavenumbers. Thus, very low Rossby numbers should be 
required in the existing proofs to guarantee that errors in the wave resonance theory are 
small for fully-developed turbulent flow. 

The present evidence from simulations and experiments is also mixed as to the validity 
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of the theory. Simulations of decaying rotating turbulence by Bardina, Ferziger & Ro- 
gallo (1985) showed a tendency for the length-scales along the axis of rotation to grow as 
rotation rate increases. Bartello, Metais & Lesieur (1994) observed two-dimensional vor- 
tices emerging from the three-dimensional flow. Hossain (1994) showed that the turbulent 
flow reduced to an approximate two-dimensional state and the energy cascaded to longer 
length scales in a 128 3 forced simulation. These results are roughly in accord with the 
wave resonance theory. The first numerical work to study explicitly the relation of flow 
two-dimensionalization and resonant triadic interactions was Smith & Waleffe (1999). 
They observed a strong two-dimensionalization and a clear inverse energy cascade. How- 
ever, they suggested that non-resonant interactions might still play an important role at 
rotation rates achieved in their simulations. So far the detailed predictions of 2D tur- 
bulence theory have not been verified numerically in rotating 3D turbulence and the 
role of resonant wave interactions has remained unclear for Reynolds number Re 3> 1. 
Thus, it is the purpose of this paper to investigate further the limit of rapidly rotating 3D 
turbulence, where resonant interactions should dominate and 2D turbulence be achieved. 

The remainder of this paper is organized as follows. In section 2, we review the resonant 
wave theory of rapidly rotating fluids, including the rigorous mathematical results. Par- 
ticular attention will be paid to a "Dynamic Taylor-Proudman theorem" that predicts 
the time-evolution of the vertically-averaged or 2D fields. In Appendix A, we give an 
elementary derivation of this result based upon the "averaged equation" of the resonant 
wave theory. In section 3, we discuss the numerical schemes and present our simulation 
results to investigate the mechanism of two-dimensionalization in rotating turbulence and 
to study the role of resonant triadic interactions. Finally, our conclusions are presented 
in section 4. 



2. Resonant Wave Theory 

Rapidly rotating fluids are a multi-time scale problem. In a rotating frame of reference, 
the Navier-Stokes equation reads (see Greenspan (1968)) 

d t u+2flxu= -VP/p+fV^-uXu. (2.1) 

Here, rotation vector $1 = £!z. u(x) is the velocity field, w = Vxu is the vorticity field, 
p is the density, v is the kinematic viscosity and P is the pressure in an inertial frame 
modified by a centrifugal term: P = P + ^p||^Xx|| 2 .IfL and U are characteristic length 
and velocity scales, then the above equation is non-dimensionalized as 

2 - - 1 - 2 

<9 t -u + — zxu = -VP + — V u - u>xu. (2.2) 
Ro Re 



Here, Rossby number Ro is defined as 



and Reynolds number Re as 



Re = (2.4) 



For simplicity, we get rid of the symbol tilde from now on. 

The first step in developing the resonant wave theory is to expand the governing fields 
into "normal modes" , i.e. the eigenmodes of the linearized equation. In the case of rapid 
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rotation this corresponds to an expansion of the velocity field into helical modes: 

u(x,t) = ^^ as (k,t)h s (k)e !k - x (2.5) 

k s=± 

where h± (k) are defined as orthogonal eigenmodes of the curl operator, satisfying ik X h s = 
s|k|h s with s = ±1 (Greenspan (1968), Waleffe (1992)). In this basis, the Coriolis term 
is diagonalized and (2.2) can be written as 

^-^ + RTe k2) ^ = \ £ EC' S '«- (2-6) 

k+p+q=0 s p ,s q 

As a convenient shorthand, we abbreviate the pair (k, s) as sk- The frequency uj Sk = 
2s(z ■ k)/fc = 2sk z /k — 2s cos 9k, with 9k the angle between Jl and wavenumber vector k. 
The zero frequency modes or so-called slow modes, with u Sk = 0, are precisely those with 
k z = and thus coincide with 2D modes having no variation along the rotation axis. The 
remaining modes with k z ^ are fully 3D modes and, since they have nonzero frequency, 
are called fast modes. The slow modes can also be obtained by vertically averaging the 
3D velocities u 3D (x, y) = jj Jq u ( x , U, z)dz, with H the vertical height of the domain. 

For Ro — ► 0, the solution of Eq.(2.6) is assumed in the resonant interaction theory 
to evolve on two distinct time-scales, the slow time t and the fast time scale r = - 
for e = Ro. A standard multiple-scale asymptotic expansion is then made with d t — ► 
d t + jfcdr- To leading order, the solution satisfies a multiple-scale Ansatz a Sk (t,r) = 
A Sk (t)e lulB k T , consisting of inertial waves with rapid oscillations on the fast time scale and 
amplitude A Sk (t) depending on the slow time t. This slow time-dependence is determined 
to eliminate secular terms growing like r, yielding the "averaged equation" : 

k+p+q=0 s p ,s q 

valid for a slow time t — 0(1). Only "resonant triads" satisfying 

^s p + u Sq + u Sk = (2.8) 

still remain in this equation. Sec Greenspan (1968), Waleffe (1993) for more details. 

There are three classes of resonant triads: "fast-slow-fast" and "slow-slow-slow" and 
"fast-fast-fast". We follow the convention that the first wavenumber in the triad is the 
one appearing on the left hand of equation (2.7), and is thus the mode which undergoes 
evolution due to interaction of the other two modes. For example, a "fast-slow-fast" triad 
gives a contribution to the evolution of a fast mode due to the interaction of another 
fast mode and a slow mode. In such a resonant triad, the slow mode acts simply as a 
catalyst for energy exchange between the two fast modes and its own energy is unchanged 
(Greenspan (1969), Waleffe (1993)). Formally speaking, there are resonant "slow-fast- 
fast" triads, but they have zero coupling coefficient. This result is known to hold more 
generally in the theory of resonant fluid interactions, not only for simple rotation but 
also for /3-plane flows (Longuct-Higgins & Gill (1967)), stratified flows (Phillips (1968), 
LeLong & Riley (1991)), rotating stratified flows (Bartello (1995)), and rotating shallow- 
water flows (Warn (1986)). As a consequence, the slow 2D modes in the limit of rapid 
rotation evolve under their own autonomous dynamics. This consists of all the "slow- 
slow-slow" triadic interactions, each of which is resonant. The averaged equation for the 
autonomous 2D modes splits into two parts, as shown by Embid & Majda (1996), Babin, 
Mahalov & Nicolaenko (1996). As Ro — > 0, the vertically-averaged horizontal velocity 
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ujf = (u 3D ,v 3D ) satisfies the 2D N-S equation: 

d t u 3 H D + (ufP .V)ul? = -VP H /p + z,V 2 u|f , (2.9) 

while the vertically-averaged vertical velocity w 3D satisfies the 2D passive scalar equation: 

d t w 3D + (u 3 H D -V)w 3D = uV 2 w 3D . (2.10) 

We give an elementary (but nonrigorous) derivation of these results in the Appendix. 
Notice that it is not implied by this result that a flow under rapid rotation will be- 
come two-dimensional, but it does mean that the dynamics will contain an independent 
two-dimensional subdynamics. In this respect, the result resembles the classic Taylor- 
Proudman theorem for steady flows (Greenspan (1968)), so that it can be termed the 
"Dynamic Taylor-Proudman Theorem." Precisely, the statement is that the "slow-slow- 
slow" triadic interactions yield the 2D-3C Navier-Stokes equations. 2D-3C means two 
variables (x,y) but three components (u 3D ,v 3D ,w 3D ). Continuing with the accounting 
of resonant triads, the final set with all fast modes are called fast-fast-fast. The averaged 
equation (2.7) for the evolution of fully 3D modes contains interactions of both the "fast- 
slow-fast" and "fast-fast-fast" types. While these cannot transfer energy directly to slow 
2D modes, Waleffe (1993) has argued that "fast-fast-fast" resonances do play an impor- 
tant role in flow two-dimcnsionalization. As a consequence of an "instability hypothesis" , 
he has shown that energy in the fast 3D modes is transferred by this set of resonances 
preferentially to slower modes, with smaller (but not zero) vertical wavenumber k z . 

The multiple-scale Ansatz and the averaged equation (2.7) have been rigorously proved 
by Embid & Majda (1996) and Majda & Embid (1998) for a general fluid dynamical 
equation with fast wave dynamics, which includes rapidly-rotating fluids as a special 
case. The precise statement of their result is that there exists some finite time T > and 
an exponent p > 1 + d/2 (with d space dimension) such that, for all < t < T, 



k,s fc 



2 



a%{t) - A.Aty*"'**'* =o(l) (2.11) 



as Ro — > 0, where a R °(t) is the solution of the full equation (2.6) for given Rossby 
number Ro, and A Sk (t) is the solution of the averaged equation (2.7). Thus the error 
in the multiple-scale approximation goes to zero in the Sobolev-norm sense of (2.11) as 
Ro — > 0. This result can be stated in another way, in terms of the error field 

S Ro (x, t) = u Ro (x, t) - U(x; t, t/Ro), (2.12) 

where u Ro (x, t) is the solution of the rotating Navier-Stokes equation (2.2) and U(x; t, t) 
is the multiple-scale Ansatz written in physical space. Define E Ro (k, t) as the wavenumber 
spectrum of the field 5 Ro (x., t), or the error energy spectrum (Kraichnan (1970), Kraich- 
nan & Leith (1972)). Then (2.11) is equivalent to the statement that, as Ro — > 0, 

/>oo 

/ dk k 2p E Ro (k,t) = o(l), (2.13) 
Jo 

for all < t < T and some p > 1 + d/2. Thus, the theorem guarantees that some moment 
of the error spectrum goes to zero, at least over a finite time-interval, as Ro — > 0. 

The multiple-scale argument applies to any rotating 3D flow, whether laminar or tur- 
bulent. However, the error bounds in (2.11) and (2.13) assume that, for d = 3, energy 
spectra decay faster than k~ 6 at high wavenumbers k, much steeper than turbulent spec- 
tra in the inertial range. Of course, when the Reynolds number Re is large but finite, 
then the spectra shall eventually decay exponentially at large enough k, in the far dissi- 
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pation range. However, because of the long inertial range, the Sobolev norms in (2.11) 
and (2.13) will become strongly Reynolds-number dependent, expected to grow as some 
power (Re)*' . Thus, the norms will not be small at high Reynolds number, except when 
the Rossby number is extremely low. Indeed, note that the spectral moments in (2.13) 
will get most of their contribution from near the Kolmogorov dissipation wavenumber 
ha in a turbulent flow. That is the last wavenumber range where resonant triads will be 
selected as Ro — > 0, because the eddy turnover time T non {kd) in the dissipation range 
is the shortest in the entire flow. Therefore, the present theorems effectively say noth- 
ing about the validity of the resonant wave theory at high Re, for realistic values of 
the Rossby number. Since previous numerical simulations of rotating turbulence by Bar- 
dina, Ferziger & Rogallo (1985), Bartello, Metais & Lesieur (1994), Mansour, Cambon 
& Speziale (1992), Hossain (1994), Smith & Walcffe (1999) have also not verified the 
predictions of (2.7), its status for turbulent flow has remained unclear. 



3. Numerical Simulations and Analysis 

To address this issue, we have simulated the rotating N-S equation Eq.(2.2) in a 128 3 
periodic box with a forcing 

/i(k,t)=e kii /ti < (k ) t)*. (3.1) 

Here, u,(k, t)* is the conjugate of Fourier component Uj(k, t). This force is specially 
chosen so that the energy input rates ek,i are all fixed (Siggia & Kerr (1978)), and we 
choose the latter to be constant in a narrow band 22 < kf < 24 and zero elsewhere. 
Thus, the total input power e = ^ k i is constant, as well as the separate energy 
inputs into slow modes and fast modes. This forcing guarantees also that the energy 
input is the same for all rotation rates. The normal viscosity term i^V 2 u is replaced 
by a hyperviscosity term (— l) p+1 v p (V) 2p u with p = 8, to extend the inertial range. 
The Coriolis force and viscosity term are integrated exactly using a slaved, 2nd-order 
Adams-Bashforth scheme. At each time step, the velocity in Fourier space u(k, t) is 
decomposed into the two helical modes and these two modes are evolved according to 
Eq.(2.6). The nonlinear term is calculated using the usual pseudospectral method and 
then projected onto the two helical modes (Smith & Waleffe (1999)). The macro-Rossby 
number Ro L = (ekj) 1 / 3 /^ ranges from 0.066 to 0.0021. Rossby numbers from the various 
experimental and numerical studies are collected in Table 1. In the table micro- Rossby 
number Ro" = u//(2Q) and u>' is r.m.s vorticity. Among the direct numerical simulations 
of forced rotating turbulence, Yeung & Zhou (1998) were interested in the dynamics in the 
range k > kf, while we and others (Smith & Waleffe (1999), Hossain (1994)) focus on the 
dynamics in the inverse energy cascade range k < kf. We consider a transient flow state, 
in which rotation is begun after a statistically steady state is reached without rotation. 
The results shown below, if without specification, were obtained after taking an ensemble 
average over eight realizations started from distinct initial conditions. This simulation 
is patterned after a 128 2 simulation of 2D N-S by Smith & Yakhot (1993). They found 
that an inertial-range energy spectrum £;~ 5 / 3 is established by the 2D inverse cascade 
process and energy flux in the inverse cascade range is a negative constant, at least 
before energy begins to accumulate at the largest scales. To compare the results of our 
3D rotating simulation with its hypothetical limit, described by the 2D N-S and passive 
scalar equations, we carry out simultaneous calculations with Eq.(2.9) and Eq.(2.10) in 
which the initial conditions are the vertically-averaged initial conditions of the rotating 
flow and the force is the vertically-averaged 3D force. In such a setup the simulation 
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Experiment/Numerical Simulations 
Traugott(1958) (decay) 
Wigeland & Nagib(l978) (decay) 
Bardina et al. (1985) (decay) 
Jacquin et al. (1990) (decay) 
Bartello et al. (1994) (decay) 
Hossain(1994) (forced) 
Yeung & Zhou(l998) (forced) 
Smith & Waleffe(l999) (forced) 
Our 128 3 DiVS'(2003) (forced) 



Ro L = 0.00064 ~ 0.0195 k f < 2 
Ro L = 0.17,0.35 fc/ = 24 



Ro L = 0.0021 ~ 0.066 22 < k f < 24 



Ro" = 1.65 None 

Ro" = 0.4 ~ 16 None 

Ro" = 0.3 ~ 6.3 None 

Ro L = 0.2 ~ 12 None 

Ro" = 0.01 ~ 100 None 



Rossby Number Forcing Scale 



Ro L = 0.1 11 < k 2 , < 13 



Table 1. Different Rossby Numbers from Experimental and Numerical Studies 



time is also the slow time. With these simulations we systematically check the main 
predictions of the averaged equation Eq.(2.7) for the resonant interactions. 

As Rossby number asymptotically approaches zero, the resonant triadic interactions 
represented in the "averaged equation" should become more and more dominant over 
the non- resonant ones. In Table 2, all resonant and non- resonant triadic interactions and 
their characteristics in the rapid-rotation limit are listed. In this section, we examine 
carefully the validity of the resonant wave theory as we decrease the Rossby number. 
To be more specific, slow-slow-slow triadic interactions, slow-fast-fast triadic interactions 
and fast-fast-fast triadic interactions are studied, respectively. We shall not study the 
catalytic "fast-slow-fast" resonant triads, because these can transfer energy only between 
fast modes with the same wavevector magnitude k and the same value of cos (Waleffe 
(1993)). Hence, they can play no direct role in transferring energy between scales or in 
two-dimensionalization of the flow. Their plausible effect is simply to isotropize the fast 
mode energy distribution in the horizontal wavenumber plane. 

In Fig. 1 are plotted 3D energy spectra for different Rossby numbers and also for the 
parallel 2D run. In the plot, kh — (k^+ky) 1 ^ 2 is horizontal wavenumber magnitude. When 
rotation is "turned on", energy is transferred to the large scales as shown in Fig. 1 (a) 
and (c), consistent with the previous observations (Hossain (1994) and Smith & Waleffe 
(1999)). Notice that large-scale energy grows faster at Ro = 0.066 than at Ro — 0.0021, 
when the energy input is the same. For Rossby number Ro = 0.17, Smith & Waleffe (1999) 
observed a rapid energy transfer to low wavenumbers similar to ours at Ro — 0.066, which 
they interpreted as due to fast, non- resonant interactions of inertial waves. For both of our 
Rossby numbers, the flow tends to two-dimensionalizc. In particular, Fig. 1 (b) and (d) 
show that slow- mode energy E(kh,k z = 0), energy from vertically-averaged horizontal 
velocity E uv (kh, k z = 0) and total energy E(k) all collapse together at large scales. This 
does not mean, however, that the fluid dynamics is that of 2D-NS. At Ro = 0.066, these 
three spectra are still far from the spectrum E2D(kh) of the 2D-NS solution. It is only 
at Ro = 0.0021 that the spectra of k z — nodes begin to approach the 2D-spectrum 
E2D(kh) suggesting that the resonant wave theory is becoming more valid for the lower 
Rossby number. 
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Type of triad Resonant/Non-resonant Characteristics as Ro — > 
slow — slow — slow R 2D + 3C 

slow — fast — fast N vanishing 

fast — slow — fast R + N catalytic 

fast — fast — fast R + N quasi — 2D 

Table 2. Triadic Interactions in Three-dimensional Rotating Turbulence 




k 

FIGURE l.Time evolution of energy spectra (upward) for (a) Ro = 0.066 and (c) Ro = 0.0021 
at t — 0, t — 172to, t = 344to, t — 430r and t — 600ro; (b) for Ro = 0.066, large-scale motion is 
strongly excited and two-dimensional over a broad range of low-wavenumbers at t = 600to (from 
one realization) and (d) for Ro — 0.0021, large-scale motion is weaker and two-dimensional over 
a smaller range at t = 600to (from one realization). E(k,to) is the initial energy spectrum. 

To investigate how energy is transferred to the large scales, we calculate the energy 
transfer functions T(kh, k z ) — ^so^) i(k z ) -Re[u*-u x cl>] at various rotation rates, where 
S(kh) denotes a circular shell of horizontal wavenumbers with central radius k^ and 
I(k z ) is an interval of vertical wavenumbers with midpoint k z . These energy transfer 
functions are plotted in Figs. 2-4, normalized by their largest absolute values. Initially 
before rotation is added, most of the energy transfer activity happens in the vicinity of 
the forcing scale fc/, where T(kh,k z ) is negative (see Fig. 2 (a) and Fig. 3 (a)). When 
the system rotates slowly, at the largest simulated Rossby number Ro — 0.066, energy 
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is quickly carried away from the forcing scales to the k z — plane and further towards 
smaller ku (see Fig. 2). There is a general tendency, at all Rossby numbers, for the 
largest transfers at later times to be into the slow modes at k z — 0. However, when the 
system rotates faster, at Ro — 0.00825, there is less energy transfer to the large scales 
at t = 450to (see Fig. 3). When the system rotates thirty times faster, at Ro = 0.0021 
(see Fig. 4), most of the energy transfer is again concentrated near the k z = plane at 
t = 450t , but hardly any transfer has developed to smaller k^. 




FIGURE 2. Normalized energy transfer functions T(kh,k z ) /\T\ max at different times for Ro = 
0.066. The red color indicates the large positive transfer and the blue color indicates the large 
negative transfer. 
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FIGURE 4. Normalized energy transfer functions T(kh, k z )/\T\ max at different times for Ro = 
0.0021. 



The resonant wave theory predicts not just that two-dimensionalization will occur, 
but also that the slow 2D modes will satisfy 2D-3C dynamics, and, in particular, should 
exhibit the phenomenology of the 2D inverse cascade. To test this, we study the spectral 
transfer at various Rossby numbers. Energy is conserved in detail for triads with all slow 
modes, so that we can define their energy flux U sss (kh) = J kh T sss (k' h )dk' h where the en- 
ergy transfer function T sss (k h ) — J2s(k h ) -Re[u*-u x <2>] only accounts for the contribution 
from slow modes k z = O.f In Fig. 5 are plotted energy fluxes H sss (kh) for different Rossby 
numbers together with the energy flux U.2D(kh) from the parallel 2D simulation. At the 
larger Rossby number Ro = 0.066, energy fluxes not only fluctuate over time but also are 
much bigger than those at smaller Rossby numbers. At Ro — 0.0021, energy fluxes are 
negative at large scales showing a complete inverse energy cascade range. Moreover, they 
begin to develop a spectral range with constant value, as expected in 2D turbulence, and 
in fact agree closely with those from 2D-NS (Fig. 5 (c)). A consistent picture is obtained 
by fitting power-laws to the energy spectra, at different rotation rates. Fig. 6 (a) shows 
slow-mode energy spectra E(kh,k z = 0) are closer to fc~ 3 than /c~ 5 / 3 at Ro — 0.066, in 
agreement with the findings of Smith & Waleffe (1999) at comparable Rossby numbers. 
However, at Ro = 0.0021 the spectrum scales as fc~ 5 / 3 (Fig. 6 (b)) similar to that of 
the parallel 2D-NS simulation in Fig. 6 (c) and as expected for a transient 2D inverse 
cascade with this forcing (Smith & Yakhot (1993)). 



t Before we calculate the transfer function T a3S (kh), we need to zero out the velocity field for 
the fast modes k z / 0. 
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FIGURE 5. Normalized energy fluxes from slow- slow- slow interactions for (a)7?o = 0.066; (b) 
Ro — 0.0021; (c) normalized energy flux from 2D Navier-Stokes at different times. Here, e 3 is 
slow-mode energy input for 3D-NS with rotation and €2n{t) is energy input for 2D-NS; ro is the 
initial large eddy turnover time. e s is the same for every rotation rate in our study. 
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FIGURE 7. Time evolution of energy spectra (upward) of (a) w 30 at Ro = 0.066; (b) w 30 at 
Ro = 0.0021 and (c) 2D passive scalar 6 2D . 

In Fig. 7 we present at various Rossby numbers the time evolution of the energy 
spectrum E w {kh, k z = 0) of the vertically-averaged vertical velocity w 3D , compared with 
the time evolution of the spectrum Eg(kh ) of a passive scalar in the 2D parallel simulation. 
Fig. 7(a) shows that at the highest simulated Rossby number Ro = 0.066 there is a 
tendency for an inverse cascade of w 3D towards the large scales. This is opposite to what 
happens for a passive scalar in a 2D inverse energy cascade range, which is well known 
to experience a direct cascade to high wavenumber: see Celani et al. (2000) and our own 
Fig. 7 (c). However, as Ro decreases, the energy spectra of w 3D become closer to those 
of the 2D passive scalar, especially at the lowest simulated Rossby number Ro = 0.0021 
(see Fig. 7 (b)). Thus, w 3D tends also to transfer downscale as the rotation rate increases. 

The difference between the 3D slow-mode dynamics and 2D-NS at finite Rossby num- 
bers is also reflected in their different flow structures. In Figs. 8-10 are plotted the 
iso-surfaces of the vertically-averaged vertical vorticity uj 3 z d from the 3D flow and those 
of the 2D vorticity w 2£) . Here, the averaged vertical vorticity is uj 3D = jj J uj z (x,y, z)dz. 
In Fig. 8 the formation of large vortices is seen at the higher Rossby number Ro — 0.066. 
Both cyclonic and anticyclonic vortices are present in the flow, which is different from 
the observations from Smith & Waleffe (1999). In their simulations, only cyclonic vor- 
tices were observed at larger Rossby numbers. At the same time, large vortices are not 
visible in the 2D turbulent flow, which seems still in its infancy (see Fig. 10). However, 
at a much smaller Rossby number Ro = 0.0021, the flow structures in 3D bear a strong 
resemblance to those seen in 2D (see Fig. 9 and Fig. 10). 

The Dynamic Taylor-Proudman Theorem predicts that the 3D slow modes will obey 
the 2D-3C Navier-Stokes equations, asymptotically for low Rossby numbers. As we dis- 
cussed in the previous section, Embid & Majda (1996) and Majda & Embid (1998) 
measured the size of the deviations from this prediction by Sobolev norms of the error 
field with large p, which do not yield useful estimates for high Reynolds number flow. 
Therefore, we shall consider here instead p — 0, or the error energy itself, as in Kraichnan 
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(1970) and Kraichnan & Leith (1972). It is useful to divide this error energy into separate 
contributions from horizontal and vertical velocity components. Precisely, we define 



As before, uff = (u 3D , v 3D ) and w 3D are the vertically-averaged horizontal and vertical 
velocities, respectively. For comparison, u 2D = (u 2D ,v 2D ) is the solution of the 2D- 
NS equation (2.9) and 9 2D is the solution of the 2D passive scalar equation (2.10). In 
Fig. 11 (a) and (b) are shown the error energies in (3.2), normalized by the energies 
En(t, Ro), Ev(t, Ro) of the corresponding 3D vertically-averaged fields, as functions of 
time at different Rossby numbers. Both plots show that the normalized error energy, 
at least over a finite interval of time, decreases as Rossby number is lowered. For any 
finite value of Ro, the vertically-averaged 3D solutions and the 2D solutions begin to 
diverge as time t increases because of their chaotic dynamics and become uncorrelated, 
leading to saturation of the normalized error energy at sufficiently long times. Although 
presently available theorems do not rigorously imply the behavior observed in Fig. 10 for 
the p = norms, nevertheless these results are quite in line with expectations from the 
formal multi-time scale analysis (Greenspan (1968), Waleffe (1993)). 



An important quantity to determine is the maximal length T* of the slow-time interval 
over which the resonant wave theory becomes valid as Ro — > . The rigorous theorems 
guarantee that a finite time T > exists so that Sobolev norms of the error (with 
p > l+d/2) converge to zero for all slow times t £ [0, T]. However, the asymptotic analysis 
does not determine the largest possible time T for which this is true. In order to estimate 
this, we have defined T(Ro) for each Rossby number Ro as the time t for which the 
normalized error energy plotted in Fig. 10(a) first reaches the value 0.5. This quantity is 
plotted as a function of Ro in Fig. 12. For a wide range of Rossby numbers above 0.00825, 
the graph can be fit by a power-law Ro^ 1 / 2 with a small error bar (only 2%). The reason 
for such a power law to exist in this range is not clear. Interestingly, there appears to 
exist a transition when Rossby number is below 0.00825 and the graph tends to saturate 
to a constant. If true, this implies that the resonant wave theory can be valid, even for 
arbitrarily small Ro, only in a finite time interval of length < T* = lim^ o ^ T(Ro). 




(3.2) 
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FIGURE 8. Time histories of the vertically-averaged vertical vorticity lJ% iso-surfaces from 3D 
rotating turbulence at Ro — 0.066. The red color is for the large positive vorticity and the blue 
one is for the large negative vorticity. 
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FIGURE 9. Time histories of the vertically-averaged vertical vorticity uj z iso-surfaces from 3D 
rotating turbulence at Ro — 0.0021. 
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FIGURE 12. Time T(Ro) vs. Rossby number Ro. There exists a power law scaling Ro 0,5 in 
the range Ro > 0.00825 but the graph appears possibly to platueau below Ro — 0.00825. 



These findings are consistent with the expected segregation of the autonomous 2D 
slow modes from the 3D fast modes at low Rossby numbers. All "slow-fast-fast" tri- 
ads transferring energy from a slow mode by interactions with two fast modes are 
non-resonant and it is another main prediction of the resonant wave theory that such 
transfer should disappear as Ro — > 0. To test this, we calculate an energy flux into 
low-wavenumber slow modes from slow-fast-fast triadic interactions. Precisely, we define 
IWfcfc) = - Io h T sff (k', k z = 0)dk' where T sff (k h , k z = 0) = T(k h , k z = 0) - T sss (k h ) 
is the energy transfer function from fast modes into slow modes at horizontal wavenum- 
ber kh- Fig. 13 (a-c) show that energy flux into slow modes at small ku decreases as Ro 
decreases. In another words, less and less fast-mode energy is transferred into the large 
scales in the 2D plane as Ro — > 0. At larger Rossby numbers, more energy is drained from 
fast modes directly into the large-scale slow modes by such non-resonant transfer, caus- 
ing the flow to become two-dimensional more quickly than at smaller Rossby numbers. 
Smith & Waleffe (1999) suggested that this non- resonant mechanism was responsible for 
the rapid two-dimcnsionalization in their simulation at larger Rossby numbers, f 



f Although we agree with their conclusion, we do not find the argument they offered very 
convincing. In their simulation, the forcing was white-noise in time with a wavenumber spectrum 
F(k) = ekexp{— 0.5(fc — kf) 2 /a 2 )/V2na 2 , nonzero in a spherical shell with mean radius k; and 
width a. They attributed the strong two-dimensionalization they observed to non-resonant inter- 
actions, because of what they claimed was purely three-dimensional forcing. However, their forc- 
ing input energy not only into the fully three-dimensional modes but also into two-dimensional 
modes in the circular band where the spherical shell intersects the k z = plane. This energy 
will inverse cascade to large scales by "slow-slow-slow" resonant interactions. 
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FIGURE 13. Normalized energy fluxes U s ff(kh) from slow-fast-fast interactions for (a) Ro = 
0.066; (b) Ro — 0.00825; (c) Ro — 0.0021 at different times. Here, e/ is the total energy input 
into the fast 3D modes. 



A final prediction of the resonant wave theory is that, as Ro — > 0, fast-mode energy 
will tend to be transferred toward smaller values of cos 9 by "fast-fast-fast" resonant in- 
teractions. Waleffe (1993) has demonstrated this tendency using the resonance condition 
(2.8) and a plausible "instability hypothesis," which asserts that energy will tend to be 
transferred out from the unstable leg of a wavevector triad. Mansour, Cambon & Speziale 
(1992) and Cambon & Jacquin (1989) have observed that energy was transferred towards 
small cos 9. However, neither distinguished fast-modes from slow-modes. Thus, the effects 
they saw may also be due to non-resonant transfer from fast modes directly into slow- 
modes. In order to investigate the "fast-fast-fast" interactions, we calculate the energy 
distribution Ef(cos9,t) only from fast modes. (We recall that the catalytic "fast-slow- 
fast" interactions can only transfer energy between fast modes with the same value of 
cos 9 and thus cannot contribute to the dynamics of this quantity.) Fig. 14 (a) shows that 
Ef(cos9,t) flattens monotonically with time at the largest Rossby number Ro = 0.066, 
implying fast modes behave largely as three-dimensional. At the middle Rossby number 
Ro = 0.00825 (Fig. 14 (b)), Ef(cos9,t) peaks at smaller cos# at earlier times and then 
flattens everywhere at later times. However, at the smallest Rossby number Ro = 0.0021, 
fast-mode energy monotonically builds up at small cos 9 (Fig. 14 (c)). This supports the 
conclusions of Waleffe (1993) based on the "instability hypothesis". 
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4. Conclusions 

To summarize briefly our results, this is the first numerical simulation of forced homo- 
geneous rotating turbulence to explore quantitatively the long-time effects of resonant 
interactions on the two-dimcnsionalization as Ro — > 0. We demonstrate an inverse cas- 

— 5/3 

cade regime with the characteristics of 2D turbulence, namely, a k h spectral range 
where energy flux is negative and constant for the slow-mode dynamics. This is con- 
sistent with the predictions of 2D-3C equations, corresponding to resonant interactions. 
Three additional findings verify the increasing importance of resonant interactions to two- 
dimensionalization at small Ro. First, vertically-averaged 3D velocities (u 3D ,v 3D ,w 3D ) 
and the solutions of 2D-3C equations converge as Ro — > 0, in the sense that the energy 
in their difference fields vanishes. Second, non-resonant energy flux into small kh in the 
k z = plane from slow-fast-fast triads decreases as Ro decreases. Finally, fast-mode 
energy is transferred toward the k z = plane, consistent with the prediction of Waleffe 
(1993) that resonant interactions of fast modes drive the flow to become quasi-2D. 

Our work has investigated rapidly-rotating turbulent fluids as a case study in the 
application of resonant wave theory. The results of the study therefore have more general 
interest than for this particular problem. The same type of theory has been widely applied 
to geophysical turbulent flows with rotation and/or stratification (Phillips (1981), Craik 
(1985)). Details of the predictions may differ considerably according to the context. For 
example, Bartello (1995) used a similar reasoning as that of Waleffe (1993) to argue 
that for rotating, stratified turbulence the catalytic, "fast-slow-fast" resonant triads will 
transfer the energy in fast gravity waves downscale and provide a primary mechanism 
of nonlinear geostrophic adjustment, while the "fast-fast-fast" resonant interactions will 
be of secondary importance. In rapidly rotating fluids without density stratification, we 
find that the opposite is true. Nevertheless, many of the basic issues are the same in 
all applications of the resonant wave theory to high Reynolds-number turbulent flows. 
Therefore, it may be useful to summarize several key conclusions from our investigation 
of this one example. 

First, the main predictions of the wave-resonance theory appear to hold for Ro <C 1 and 
Re ^> 1 simultaneously. Convergence of the approximation holds in the useful sense that 
the energy in the error fields decreases to zero, over a finite time-interval. Present rigorous 
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theorems (Embid & Majda (1996), Embid & Majda (1998), Majda & Embid (1998)) use 
Sobolev-norm estimates that yield useful estimates for rotating laminar flows, but not for 
fully-developed turbulence. An optimistic conclusion of our work is that those theorems 
can probably be improved, by substitution of an L 2 or energy norm for the Sobolev 
norms. A proof of this may not be easy, given the current lack of understanding of the 
inviscid, high-Reynolds number limit of fluid equations. However, our numerical results 
lend support to the idea that the energy-containing, large-scale motions in turbulent flow 
will conform to the predictions of resonant wave theory in the appropriate limit. 

We find that there is an intrinsic time restriction on the validity of the resonant wave 
theory. The formal multiple-scale analysis, as well as the rigorous theorems, predict that 
there should be a finite interval of the slow-time variable t over which the approximation 
will converge. They do not answer the question whether this time-interval is finite in 
fact, or whether the approximation will be valid at any arbitrary time for a sufficiently 
small Rossby number. Our numerical results seem to indicate either that the convergence 
holds over only a finite time-interval or else that the Rossby numbers required for a 
close approximation become extremely small beyond a certain time-horizon. Of course, 
both the dynamics involved — rotating 3D-NS (2.2) and the low Rossby-number limiting 
dynamics given by the averaged equation (2.7) — are chaotic at high Reynolds number and 
a rapid exponential divergence of their solutions is to be expected for finite Ro. This does 
not imply a necessary failure of the resonant wave theory applied to forced, statistically 
stationary turbulence, since chaotic dynamics often possess "structural stability" of their 
long-time statistics. 

A more serious limitation of the wave-resonance theory may be the magnitude of 
the Rossby numbers required for convergence, even at relatively early times. We have 
found sizable non-resonant effects for Rossby numbers as small as Ro = 0.066. Since 
the Rossby number is not extremely small in most engineering applications, flow two- 
dimensionalization is there probably due mainly to the non-resonant interactions, for 
which the 2D dynamics is not segregated from the 3D dynamics. However, the resonant 
wave theory provides the correct "rapid-distortion limit" for fast rotations. Therefore, it 
may still be useful to take into account the resonance conditions in statistical modeling 
schemes for engineering purposes, as a constraint which becomes exact as i?o ^ 0. 
Resonant interactions may be more important for geophysical flows, where typical Rossby 
numbers are much smaller. However, for the atmosphere at midlatitude synoptic-scales, 
the typical Rossby number is only about Ro = 0.1 and even for the Gulf Stream with 
smaller characteristic velocities the Rossby number is only Ro = 0.07 (Pedlosky (1987)). 
In our simulations, this is a marginal value to observe the predictions of the theory A 
combination of rotation with other effects such as stratification may be necessary in the 
atmosphere and oceans to select resonant interactions of fast waves. 
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Appendix A. Dynamic Taylor-Proudman Theorem 

Here, we use the helical decomposition to give a simpler and self-contained derivation 
of the Dynamic Taylor-Proudman Theorem for the "slow-slow-slow" resonant triadic 
interactions. In a triad consisting of all slow modes k,p and q, then k z — 0, p z — and 
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q z = 0. In the k z = plane or 2D plane, the wave number vector k = k x x + k y y and 



its amplitude k = + k y . By choosing the helical modes h s (k) = Nxk + isv which 
satisfy ikxh± = ±|k|h± (Waleffe (1992)), we have 

h s (k) =z + is ky± ~ kxf =z + z,sk ± . (Al) 

Here N = and k- 1 = (fc y x — k x y)/k are unit vectors orthogonal to the unit 

vector k = k/fc. The other two helical modes in the triad are h s (p) = z + isp 1 - and 
h s (q) = z + isq 1 - , where the unit vectors p x = (p y x — p x y/p and q^ = (q y x — q x y)/q 
are perpendicular to the unit wave vectors p = p/p and q = q/g, respectively. The 
velocity projection onto h s is 

a (k k t) — ^sfe' 0) ' "(^" 2j 

= l[u z (k x ,k y ,0)+is(^-^)}. (A2) 

Since 

w = ikxu = i(fca;X + k y y) X (u^x + u y y + u z z) = i(k x u y — k y u x )z = t£> z z, (A 3) 
Eq.(A2) can be also written as 

a s {k x ,k y ,0,t) = ^u z {kx,k y ,0) + f ^feA' ) . (A4 ) 

Thus, the vertically-averaged vertical velocity and the vorticity can be separated from 
the above equation as 

ii z = a + + a-, uj z = k(a + — a_). (A 5) 

For the slow modes a Sk (t) = A Sk (t), the "averaged equation" Eq.(2.7) becomes 



(d t + ± k 2 )a Sk = \ E ( s pP s q q)(K p xK t ) ■ K k a; p a* Sq . (A 6) 

Since ^Xq 1 = p X q, 



Re ' Sk 4 

k+p+q— s p ,s q 



(K p XK q )-K k = (z-«sp ± )x(z-«sq ± ) • (z-isk 1 ) 

= (s p p - s fc k) X (sfck - s 9 q) • z . (A 7) 

Then, Eq.(A6) becomes 

(dt + ^ 2 )a Sfc = j E E( s P p ~ s ^)(%>P- s fck)x(s fe k-s 9 q)-z. ( A8 ) 

k+p+q=0 Sp.s, 

To construct u z and u) 2 (Eq.(A5)) from the above equation, we have 

( d t + k 2 )(a +k - a_ k ) = ^ ^ ^_^_( pX q • z)[-a* +p a* +q - a*_ p a*_ q 

k+p+q=0 



and 



+a* +p a*_ q +a*_ p a* +q ], (A 9) 

(dt + ^;k 2 )(a +k +a- k )=\ £ (pxq • z)[(p - q)(~a* +p a* +q + a*_ p a*_ q ) 

+ (p + q)K p a*_ q -a*_ p a* +q )} .(A 10) 



Re V ~rfc — 2 

k+p+q=0 



22 Q. Chen, S. Chen, G. L. Eyink and D. D. Holm 

Let b k = a +k -a_ fc and 9 k = a +k + a_ fc , then a +k = (6 k + b k )/2 and a_ fc = (6 k -b k )/2. 
Substituting the above relations into Eq.(A9) and Eq.(AlO), we have 

{dt + h k2)bk = ~\ £ (ah) 

k+p+q=0 

and 

(9 t + ^m = i £ ^^W^W ^ 12) 

k+p+q=0 Fq 

Since t2> z = kb k from Eq.(A 5), Eq.(A 11) can be written in terms of the vorticity as 



(3 t + -fc 2 )^- £ i^_[^- P 6^;]. (A 13) 

k+p+q=0 ^ 

This is exactly the same as the helical formulation of 2D N-S equation which is given by 
Waleffe (1993). 

Also because of u z (k) = 9 k from Eq.(A5), Eq.(A 12) can be written in terms of u z (k) 

as 

k+p+q=0 ^ 

This equation has the same structure as the 2D vorticity equation Eq.(A 13). It is known 
that a 2D passive scalar and the 2D vorticity obey homologous equations and, in fact, 
the above equation is the helical formulation of the 2D passive scalar equation. 

Based on the above discussion, we now summarize the Dynamic Taylor-Proudman 
theorem as the following: in three-dimensional rotating turbulence, when Ro — > 0, "slow- 
slow-slow" triadic interactions yield the 2D-3C Navier-Stokes equations. Here, 2D-3C 
refers to two variables (x,y) and three components (u 3D ,v 3D ,w 3D ). The vertically- 
averaged horizontal velocity (u 3D ,v 3D ) = jj J^ 1 (u x (x, y, z),u y (x, y, z)dz satisfies the 2D 
Navier-Stokes equation 

d t Lb z + vk 2 uj z = ~ik x lI 3D (x,y)uj z (x,y) - ik y v 3D (x, y)u) z (x, y), (A 15) 

and the vertically-averaged vertical velocity w 3D = u z (x,y, z)dz obeys the 2D 

passive scalar equation 

dtW^ik) + vk 2 W 3b (k) = -ik x u 3D (x,y^ 3D (x,y)-ik v v 3D (x,yyw 3D (x,y). (A 16) 

Here, H is the vertical height of the domain and the three averaged velocity components 
in the Fourier space are as follows: 

^( k )= i TCT^( k )'^( k ) = - i 7CT^( k )' ^( k )=^( k )- (A17) 

K x "r K y x K y 
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